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SUMMARY 

This  document  describes  the  theory  and  input  requirements  to  the 
CRKGRW  Fortran  computer  program,  which  is  used  to  calculate  crack  growth 
rates  using  data  on  material  properties,  load  sequences  and  crack  geometry. 
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NOTATION 


a  -  crack  length  or  depth  (see  Figs  4,  5  and  6) 
a, i  -  initial  crack  length  or  depth 
B  -  constant  in  modified  Forman  equation 
b  -  width  or  half  width  of  specimen  (see  Figs  4,  5  and  6) 
C  -  constant  in  crack  growth  equation 
Cp  -  Wheeler  retardation  factor 
c  -  crack  depth  or  length  (see  Figs  4,  5  and  6) 

Co  -  initial  crack  depth  or  length 
da/dN  -  constant  amplitude  crack  growth  rate 
dc/dN  -  constant  amplitude  crack  growth  rate 
K  -  stress  intensity 
K,  -  fracture  toughness 
Kjn  -  threshold  stress  intensity 
A  A'  -  change  in  stress  intensity 
L  -  characteristic  length 
M  -  constant  in  crack  growth  equation 
m  -  Wheeler’s  shaping  exponent 
n  -  exponent  in  crack  growth  equation 
t*  -  constant  in  modified  Forman  equation 

Q  -  constant  in  modified  Forman  equation 

K  -  stress  ratio  ie  c,,,,,. /a,„„x 
Rfl  -  extent  of  current  yield  zone 

r  -  hole  radius 

t  -  thickness  or  half  thickness  of  specimen 
(see  Figs  4,  5  and  6) 

P  ~  geometrical  factor 

a  -  current  stress  level 

cr„/,  -  cycle  alternating  stress 

o max  ~  cycle  maximum  stress 
a,„e„, ,  -  cycle  mean  stress 

cr,„j„  -  cycle  minimum  stress 
crrfl)  -  reduction  in  stress  for  Willenborg  model 
<jy  -  yield  stress 

Ait  -  cyclic  stress  range 


DEFINITIONS 


The  terms  LAYER,  MISSION,  FLIGHT,  SEGMENT  and  BLOCK  are  used  in  this 
manual.  Their  definitions,  as  they  apply  in  this  manual,  are  set  out  below  (also 
refer  to  Fig.  8). 


LAYER  -  A  layer  consists  of  one  or  more  constant  amplitude  loading 
cycles. 

MISSION  -  A  mission  is  a  group  of  layers  defined  in  a  specified  manner. 

FLIGHT  A  night  is  oxic  application  of  a  mission. 

SEGMENT  -  A  segment  consists  of  one  or  more  flights  of  the  same  mission 
applied  consecutively.  The  flights  within  a  segment  are  num¬ 
bered  from  one  upwards. 

BLOCK  -  A  block  consists  of  all  the  segments  specified  in  the  data  file. 


Note:  Random  sequences  can  be  handled  by  CRKGRW  if  each  cycle 
is  considered  as  a  separate  layer. 
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1.  INTRODUCTION 


This  report  is  a  manual  for  using  CRKGRW  Version  2.0  (hereafter  called  CRK- 
GRW)  to  calculate  crack  growth  using  various  models.  Relevant  background  the¬ 
ory  is  also  included. 

CRKGRW  is  essentially  an  upgrading  of  the  FORTRAN  77  program  CRKGRW 
Version  1.0  with  program  bugs  removed  and  more  options  added.  CRKGRW 
Version  i.O  is  itself  a  modification  of  CRACKS  IV.  This  was  a  FORTRAN  IV 
program  written  by  R.  M.  Engle  Jr  of  the  Air  Force  Flight  Dynamics  Laboratory, 
Wright  Patterson  Air  Force  Base. 

This  manual  also  includes  elements  from  two  ARL  reports  on  CRACKS  IV  [Refs 
11  and  12). 

2.  CRACK  GROWTH  MODELLING 

The  essence  of  crack  growth  programs  is  to  calculate  the  change  in  crack  length 
'■esulting  from  the  application  of  a  single  load  cycle.  This  is  repeated  for  each  load 
cycle  in  the  sequence  resulting  in  the  total  crack  length  and  the  time  to  failure 
being  calculated.  / 

The  basic  material  data  used  for  the  calculations  are  generally  available  in  the 
form  of  a  da/dN  versus  A K  relationship  (see  Section  3.2).  In  order  to  calculate 
the  increment  of  crack  growth  per  load  cycle,  the  value  of  the  change  in  the  stress 
intensity  factor  A K  at  the  crack  tip  must  be  determined.  This  factor  requires  a 
knowledge  of  the  crack  length,  specimen  and  crack  geometry,  and  the  magnitude 
of  the  load  cycle  so  that  a  /?  factor  can  be  calculated.  For  linear-elastic  fracture 
mechanics  the  equation  relating  A K  with  the  above  factors  is  given  in  Equation 
la. 


A/f  =  Aas/iva(l  (la) 

This  equation  does  not  take  into  account  the  subsequent  effect  of  large  stress  cycles 
which  may  retard  the  growth  rate.  This  can  be  incorporated  in  a  number  of  ways 
(as  described  in  Section  3.3). 

CRKGRW  calculates  crack  growth  layer-by-layer  according  to  the  equation  below. 

a,  —  a,_ i  +  {da/dN),  *  (Number  of  cycles  in  the  layer)  (16) 

where  [da/dN),  is  crack  growth  rate  for  a  single  cycle  in  the  layer. 

If  failure  occurs  within  the  layer,  then  CRKGRW  will  re-calculate  the  crack  growth 
in  the  layer  cycle-by-cycle.  In  this  case  only,  the  crack  length  is  updated  after 
each  cycle  until  failure  occurs. 


3.  THEORETICAL  BACKGROUND  TO  PROGRAM 


3.1  Introduction 

This  section  has  been  included  so  that  the  user  of  the  CRKGRW  program  will  have 
some  theoretical  knowledge  as  to  how  the  program  calculates  various  parameters. 
Crack  length  a  and  depth  c  for  each  solution  will  be  defined  in  Section  3.4. 

3.2  Material  Data 

The  material  data  required  by  CRKGRW  are: 

i)  material  yield  stress, 

ii)  fracture  toughness, 

iii)  stress  intensity  threshold,  and 

iv)  da/dN  versus  A K. 

The  first  three  are  simply  numbers  but  the  fourth  can  take  many  forms.  CRKGRW 
allows  several  different  da/dN  versus  A K  relationships  and  these  are  described  in 
the  next  section. 

3.2.1  da/dN  versus  A K 

Five  equations  are  presented  which  attempt  to  model  the  behaviour  of  the  da/dN 
versus  AK  curve  (Fig.  1)  by  taking  various  parameters  into  account. 

a)  Paris  equation  (Ref.  1] 

This  equation  has  the  following  form:- 

da/dN  =  C{AK)n  (2) 

The  program  requires  the  specification  of  the  parameters  C  and  n. 

b)  Forman  equation  [Ref.  2] 

The  Forman  equation  has  the  following  form:- 

C(AK)" 

dN  (1  -  R)Kr  -  AK  (  > 

This  equation  takes  into  account  the  influence  of  the  of  the  stress  ratio 
R  and  the  material  fracture  toughness  Kc. 

c)  Forman  equation  modified  (Ref.  3) 

The  modified  Forman  equation  has  the  following  form:- 

da  _  C{AK  -  [PR  +  Q))AK" 
dN  ~  (1  -  BR)Kr  -  AK 

This  equation  has  been  demonstrated  to  reproduce  Fig.  1,  for  different 
stress  ratios,  very  accurately.  The  5  parameters  required  are  C,  n,  P,  Q 
and  B. 

d)  Walker  equation  (Ref.  4| 

Walker’s  equation  in  original  form  was  an  extension  of  the  Paris  equation 
to  include  the  influence  of  the  stress  ratio  R.  CRKGRW  uses  the 
following  form:- 


J 


The  three  parameters  C,  M  and  n  are  required  by  CRKGRW. 
e)  Tabulated  data 

Often  the  most  convenient  form  of  input  of  constant  amplitude  data  for 
da/dN  versus  A K  is  tabulated  values.  These  are  obtained  from  graphs 
such  as  Fig.  1. 

3.2.2  Evaluation  of  Material  Constants 

Representative  material  data  having  been  selected,  (either  from  handbooks  or 
from  constant  amplitude  test  results),  they  now  need  to  be  fitted  to  the  equations 
presented  in  section  3.2.1.  The  Paris  and  Forman  constants  C  and  n  are  sometimes 
quoted  in  materials  data  books.  This  section  looks  at  the  determination  of  the 
modified  Forman  equation  constants  as  an  example. 

i)  Evaluation  of  P  and  O.  P  and  Q  describe  the  threshold  stress  inten¬ 
sity  factor  range  and  its  movement  with  R.  Figure  2(a)  shows  how  the 
constants  can  be  obtained  from  the  da/dN  data,  if  more  than  two  stress 
ratios  are  given. 

ii)  Evaluation  of  B.  8  describes  the  variation  of  the  fracture  toughness 
with  R.  Figure  2(b)  shows  how  B  can  be  evaluated  from  the  da/dN 
data. 

iii)  Evaluation  of  C  and  n.  C  and  n  are  obtained  from  the  line  of  best  fit 
passing  through  the  approximately  linear  centre  section.  C  is  calculated 
from  the  vertical  intercept  and  n  from  the  slope.  Some  compromise  may 
be  necessary  when  evaluating  C  and  n  to  obtain  the  best  fit  over  the  full 
range  of  R. 


da  _  r  A  K 

dN  ~  '  [(1  -  fl)t‘--''7 


3 


3.3  Retardation  Models 

Retardation  models  take  into  account  the  reduction  in  crack  growth  rate  after 
the  application  of  a  high  load.  At  present  the  CRKGRW  program  contains  two 
working  retardation  models.  These  models  are:- 

a)  Wheeler  [Ref.  5] 

This  model  introduces  a  retardation  factor  Cp  which  reduces  the  growth 
rate  function,  /,  as  described  by  the  relationship 

n 

a„  =  aft  -f  Y,Cr.f(*K,)  (6) 

where 

cr.  =  for  (a+  R„<  ar) 

or 

Cpt  =  1  for  (a  +  Ry  >  ap)  ie:  no  retardation 

and 

ap  -  a  =  the  distance  from  the  current  crack  tip  to  the  furthest  point 
of  the  elastic-plastic  interface. 

The  parameter  m  must  be  determined  experimentally  or  obtained  from 
data  handbooks.  The  parameter  is  also  a  useful  tool  in  calibrating  the 
program  run  to  test  data. 

b)  Willenborg  [Ref.  6] 

Willenborg  modelled  retardation  by  reducing  both  the  maximum  and 
minimum  stress  of  cycles  following  an  overload.  The  reduction  in  these 
stresses  is  given  by> 

_  _  f./2(ap-a)  _  ,  , 

Or,  , I  ~  p  y  a  Otnrtr  t  ‘) 

provided  that  <T,„al  has  been  preceded  by  a  higher  load,  <7!  which  pro¬ 
duced  ap 

This  retardation  model  has  the  effect  of: 

i)  reducing  the  stress  ratio  R  leaving  A K  unchanged;  or 

ii)  reducing  A/f  for  R  =  0;  or 

iii)  making  AA'  and  R  both  equal  zero  when  (Ti/a,,,,,,  >  2.0 

In  its  standard  form  this  model  requires  no  additional  input  parameters. 
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3.4  Stress  Intensity  Solutions 

Every  crack  growth  program  is  required  to  calculate  the  stress  intensity  solution. 
CRKGRW  has  the  capability  to  handle  several  different  solutions  for  fi  .  These 
solutions  are  described  below. 

Note:  In  the  solutions  below,  if  no  definition  of  a  and/or  c  is  given  then: 
i)  the  problem  is  regarded  as  one  dimensional  and 
ii)  a  is  defined  as  along  the  surface  of  the  specimen. 

a)  Surface  Flaw  Solution 

The  surface  flaw  solution  written  into  CRKGRW  is  the  Newman  and 
Raju  solution  | Ref.  7,  Equation  17].  Figure  4  shows  the  definition  of  a 
and  c.  The  equation  for  p  is  given  below:- 

P  =  Mt  +  j  +Mj^-^  <;/o/u  (**) 

provided  the  following  conditions  ^re  satisfied 
0  <  c/a  <  2  a/b  <  0.5,  0  <  </>  <  jr 

where  <p  =  0  for  a  and  <*>  =  |  for  c. 

When  the  crack  becomes  a  throngh-crack  ft  is  automatically  made  unity 
since  the  problem  is  now  a  through-crack  in  a  plate. 

Definition  of  the  variables  can  be  found  in  Ref.  7,  pages  12-13. 

Note:  -  This  solution  includes  its  own  finite  width  correction  factor. 

-  a  in  this  manual  is  the  same  as  c  in  Ref.  7  and  c  in  this 
manual  is  the  same  as  a  in  Ref.  7. 

b)  Constant  value  of  ft 

BETA  option  1  allows  the  user  to  input  a  constant  value  of  ft  into  the 
program. 

c)  Finite  width  correction 

BETA  option  2  allows  for  a  finite  width  correction.  The  correction  used 
was  developed  by  Fedderson  jRei.  sj.  The  equation  U  shown  below 


where  b  is  the  half  width  of  the  plate  and  r  is  the  radius  of  the  hole,  if 
any,  in  the  plate  . 

This  correction  may  be  used  with  any  of  the  other  solutions  except  for 
the  Surface  Flaw  solution  and  BETA  options  10  and  11. 

d)  Tabular  input  of  ft  solution 

BETA  options  3  and  4  allow  tables  of  P  versus  <2/ L  to  be  entered. 

e)  Bowie  solution  for  a  single  crack 

BETA  option  5  is  the  Bowie  solution  for  a  single  crack  from  a  circular 
hole  [Ref.  14,  Table  lj.  The  equation  is  shown  below:- 


5 


0.8734 


(10) 


P 


0.0702  + 


(0.3240  +  «) 


f)  Bowie  solution  for  a  double  crack 

BETA  option  0  is  the  Bowie  solution  for  a  double  crack  from  a  circular 
hole  (Ref.  14,  Table  lj.  The  equation  is  shown  below:- 


P  =  0.9439  + 


0.0805 
(0.2772  +  4) 


HI) 


g)  ASTM  compact  tension  specimen 

BETA  option  7  allows  the  solution  for  ASTM  compact  tension  specimens 
(Ref.  10,  Equation  12].  This  equation  is  shown  below  (see  also  Fig  7). 
Note  that  the  spectrum  is  input  as  loads  and  not  stresses. 


where 


P  = 


Jil.)= 

t\fiab 


/  (!)  .  4.55  .  <0.32  (!)  K,4.,  (!)'  -  IW»  (!)’  7  37»,  Q' 

-  4287(2)%  2017  (!)‘ 
for  0.2  <  a/W  <  0.8 

h)  Grumman  compact  tension  specimen 

BETA  option  8  allows  the  Grumman  compact  tension  specimen  to  be 
modelled  [Ref  13,  Equation  40].  See  Fig.  7  for  the  geometry  of  the 
specimen.  The  equation  for  /?  is: 

(0.1229  +  16.4098(£)  ~  37.395( )2  +  54.7007 (£)3) 

P 

for  a/b  <  0.5 


(114.054  -  830.132(2)  +  2327.177 (4)2  -  2890.81 1  ( if  +  1382. 200(  £)4) 

s/iTa/t 


for  a/b  >  0.5 

i)  Liu  crack  from  a  circular  hole 

BETA  option  9  allows  the  Liu  solution  for  a  double  quarter  crack  from 
a  circular  hole.  The  equation  used  is  shown  below:- 


1.12  -  -  0.0702  + 


0.8734 
(0.3246  I- 
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Note  the  correspondence  between  this  solution  and  the  Bowie  solution  for  a 
single  crack. 

j)  Corner  crack  from  a  hole 

BETA  option  10  allows  th»  Newman  and  Raju  solution  for  a  corner  crack 
from  a  hole  (Ref.  7].  Figure  5  shows  the  definition  of  a  and  c.  The 
equation  for  [3  is  given  below:- 

q  4 1 

p—  +  +  Mi  I  5i 9293  fofw  (13) 

provided  the  following  conditions  are  satisfied 

0.2  <  c/a  <  2  0.5  <  c/t  <  1,  (r  -t  a)/b  <  0.5,  0  <  <fi  <  | 

where  ij>  =  0  for  a  and  <f>  =  |  for  c. 

Definition  of  the  variables  can  be  found  in  Ref.  7,  pages  20-22. 

See  also  the  note  for  the  Surface  flaw  solution. 

k)  Surface  crack  in  a  hole 

BETA  option  11  allows  the  Newman  and  Raju  solution  for  a  single  or 
double  embedded  crack  in  a  hole  [Ref.  7].  Figure  6  shows  the  definition 
of  a  and  c.  The  equation  for  0  is  given  below:- 

P  =  Mi  +M2(^)  j  gig?gif0f\y  (H) 

provided  that  the  same  conditions  given  for  BETA  option  10  are  satisfied. 

Definition  of  the  variables  can  be  found  in  Ref.  7,  pages  10-20. 

See  also  the  note  for  the  Surface  flaw  solution. 


4.  DATA  INPUT  TO  CRKGRW 


The  input  data  file  contains  a  number  of  headings,  called  KEYWORDS,  followed 
by  the  relevant  information  on  the  next  line  or  lines  for  that  keyword.  This  section 
describes  some  areas  wheie  care  must  be  taken  to  produce  an  error  free  data  file 
and  then  a  description  of  all  the  possible  data  inputs  follows. 

UNITS 

Any  consistent  set  of  units  may  be  used. 

COMMENT  CARDS 

A  C  at  the  beginning  of  a  line  represents  a  comment  card.  Care  must  be 
taken  not  to  place  comment  cards  where  data  is  expected  by  the  program. 
Comment  cards  are  only  allowed  before  a  keyword  or  after  all  information 
associated  with  the  current  keyword  has  been  input. 

DEFAULT  VALUES 

For  some  of  the  options  discussed  in  the  following  section,  if  a  zero  is 
specified  then  a  default  value  will  be  used  (which  may  be  non-zero). 

KEYWORDS 

CAPITALS  should  always  be  used  in  the  data  file  fcr  the  keywords.  The 
program  will  not  recognise  lower  case.  Lower  case  may  be  used  when 
writing  a  comment  or  heading  line  only. 

All  keywords  should  be  placed  at  the  beginning  of  the  line  (ie:  they  must 
begin  in  column  l). 

ORDER  OF  KEYWORDS 

Some  keywords  may  be  placed  anywhere  in  the  data  file  while  others 
must  occur  in  a  certain  order.  To  avoid  problems,  use  the  order  which  is 
specified  in  section  4.1  . 

OPTIONAL  KEYWORDS 

Not  all  of  the  keywords  need  to  appear  in  the  data  file.  Those  keywords 
that  are  optional  are  indicated  in  section  4.1  as  they  are  discussed. 

NUMERICAL  INPUT 

All  the  numerical  values  required  by  the  keywords  may  be  entered  in  free 
format.  That  is,  the  values  may  appear  anywhere  on  the  line  in  the  order 
that  they  are  required  with  either  a  comma  or  a  space  separating  them. 

ERRORS  IN  INPUT 

The  program  has  been  designed  to  recognize  some  of  the  errors  that 
may  oc'-ur  in  the  input  file.  However,  it.  will  not  catch  sillv  errors  such 
as  negative  dimensions.  Some  of  the  m<  f  common  errors  are  given  in 
section  4.'(. 


T" 


} 


4.1  KEYWORDS 

The  keywords  recognised  by  CRKGRW  ns  flagging  a  data  segment  are: 
LINES  PER  PAGE 
TITLE 
EQUATION 
MATERIAL 
THRESHOLD 
LIMITS 
ANALYSIS 
LOADS 
SPECTRUM 
PRINT 
PLOT 
END  DATA 


4.1.1  LINES  PER  PAGE  segment  (optional) 

This  segment  defines  the  number  of  lines  of  output  that  will  be  printed  per  page. 
The  default  is  80. 

Line  1:  LINES  PER  PAGE 
Line  2:  L 

L  =  the  number  of  lines  per  page  (must  be  an  integer  >  40) 


4.1.2  TITLE  segment  (optional  but  should  be  used) 


Line 

1: 

TITLE 

Line 

2: 

n 

where  n  is  the  number  of 

Line 

3: 

p, 

line  of  title  information. 

Line 

4: 

2n<! 

line  of  title  information. 

Line 

n  *4-  2 1 

n'1' 

line  of  title  information. 

The  title  text  will  be  printed  at  the  head  of  the  output  file.  A  maximum 
of  80  characters  is  allowed  per  line. 
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4.1.3  EQUATION  segment 

This  segment  determines  the  way  the  da/dN  versus  A K  data  will  be  input  to  the 
program.  There  are  eight  different  equations  the  program  will  accept.  They  are:- 
PARIS 

PARIS-B1  (bi-linear  version  of  the  Paris  equation) 

FORMAN 

FORMANMOD 

WALKER 

WALKER-BI  (bi-slope  version  of  the  Walker  equation) 

DA/DN 

R-DA/DN 

CRKGRW  allows  two  definitions  of  stress  intensity,  K  :• 

i)  K  =  ffy/rraf} 
and 

ii)  K  =  c\/afd  (NASA  convention) 

CRKGRW  will  default  to  the  first  definition  unless  the  keyword  NASA  is  placed 
after  the  above  equation  options. 

Line  1:  EQUATION 

Line  2:  Characters  1  to  10  -  one  of  the  above  equation  keywords. 

Characters  11  to  14  -  NASA  (if  required). 

4.1.4  MATERIAL  segment 

The  nature  of  input  in  this  segment  is  determined  by  the  option  previously  chosen 
in  the  EQUATION  segment. 

Line  1:  MATERIAL 

Line  2:  This  is  a  comment  line,  which  is  used  to  describe  the  material  in  use 
(maximum  of  80  characters). 

It  is  written  to  the  output  file. 

Line  3:  This  line  provides  the  data  for  the  EQUATION  chosen. 

a)  PARIS  option 

Line  3:  C  and  n 

b)  PARIS-BI  option 

Line  3:  C \  ,  nj  ,  A KtRaxs  i  U*2  j  n? 

This  gives  a  bi-linear  Paris  fit  with  the  parameters 
C\  and  ri|  for  A K  <  AKrn.\y>' 
and  On  and  ni  for  A K  >  A  Kti;a  xs  ■ 

c)  FORMAN  option 

Line  3:  C,  n,  K, 

d)  FORMANMOD  option 

Line  3:  C,  n,  K,-,  P,  Q,  B 

e)  WALKER  option 
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Line  3:  C,  M,  n 


f)  WALKER-BI  option 

Line  3:  Ci,  Mi,  nj 
Line  4:  <72,  A/2,  ^2 
Line  5:  ( da/dN)ri!A\s ,  A KtjiaXS 

Figure  3(a)  shows  the  set  of  lines  generated  by  the  standard  Walker 
equation  for  a  range  of  R  values.  Note  that  there  are  no  lines  for  R  <  0 
as,  at  present,  CRKGRW  considers  all  negative  loads  to  be  zero. 
Figure  3(b)  shows  the  effect  of  the  WALKER-BI  option.  As  can  be 
seen  this  splits  the  Walker  equation  into  two  regions  such  that: 

Ci,  Mi,  and  ni  apply  in  region  I 
and  C2,  A/2,  and  ri2  apply  in  region  2 

where  the  transition  between  the  two  regions  is  defined  by: 
da/dN  =  (da/ dN)TR.\\ s  and  A K  =  A Ktpa.xs  on  the  R  =  0  line. 

g)  DA/DN  option 

A  table  of  values  of  AK  versus  da/dN 

Line  3:  n  -  the  number  of  points  in  the  table  (must  be  an  integer) 
Line  4:  l'1  point,  A Ki,  ( da/dN)\ 

Line  5:  2"1'  point,  AK2,  (da/dN) 2 

Line  n  +  3:  n'1'  point,  A K„,  (da/dN), , 

If  line  3  is  negative  then  CRKGRW  assumes  that  lines  4  to  n+3  are 
in  a  separate  data  file.  The  name  of  this  data  file  will  be  requested  by 
CRKGRW  when  it  is  needed. 

h)  R-DA/DN  option 

Sets  of  constant  amplitude  crack  propagation  data  for  up  to  ten  values 
of  stress  ratio  R  may  be  entered.  The  data  for  each  value  of  R  must  be 
placed  in  separate  files.  File  names  will  be  requested  at  runtime  and 
must  be  given  in  order  of  increasing  R.  'I'll"  format  of  each  file  is  as 
follows:- 

Line  1:  R  value 

Line  2:  n  -  the  number  of  points  in  the  table  (must  be  an  integer) 
Line  3:  1st  point,  A/f),  (da/dN)  1 
Line  4:  2"'1  point,  A  Ki,  (da/dN)2 

Line  n  +  2:  n,K  point,  A K„,  (da/dN)„ 

The  final  line  of  the  MATERIAL  segment  requires  the  parameters  K,, 
and  a,,. 

Kq  -  maximum  permissible  value  of  stress  intensity. 
a ,,  -  material  yield  stress. 
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4.1.5  THRESHOLD  segment 

Line  1:  THRESHOLD 
Line  2:  AK-j-h,  R.viilt 

AKjh  ~  is  threshold  value  of  A K. 

Rmiilt  ~  parameter  if  crack  growth  is  R-dependent 
Default  =  0 

Crack  increments  will  be  calculated  if 

K  >  AKth{1  -  *  R) 

4.1.6  LIMITS  segment 

Line  1:  LIMITS 

Line  2:  a„  ,  am„,  ,  An  ,  Rout 

a, i  -  Initial  length  or  half  length  of  crack.  Care  must  be  taken  with  the 
definition  of  an.  Refer  to  the  theory  section  to  determine  the  definition 
of  a  pertaining  to  the  /?  solution  being  used. 

a, i, nr  -  Maximum  length  or  half  length  of  crack. 

Default  =  102" 

N„  -  Initial  cycle  number.  Cycles  applied  during  the  run  will  be  numbered 
starting  with  this  value. 

Default  -  1 

Rear  ~  Calculated  stress  ratios  are  truncated  to  this  value. 

Default  =  1, 
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4.1.7  ANALYSIS  segment 

Several  options  are  available  for  analysing  the  stress  intensity  for  various  specimen 
geometries.  Two  or  more  of  these  options  may  be  selected  to  operate  simultane¬ 
ously  and  different  options  can  operate  over  different  crack  lengths. 

Line  I:  ANALYSIS 

An  option  is  selected  by  choosing  one  of  the  keywords  below,  followed  on  the  next 
line  by  the  various  numerical  parameters.  Valid  keywords  are:- 

SURFACE 

RETARD 

BETA 

END 

One  or  more  of  the  keywords  SURFACE,  RETARD  and  BETA  must  be  used. 

a)  SURFACE  option 

Line  1:  SURFACE 
Line  2:  Cl  ,  C2  ,  C3  ,  C4 

Cl  -  initial  crack  depth,  c. 

Default  =  a„. 

C2  -  thickness  of  plate,  t. 

C3  -  width  of  plate,  2b. 

C4  -  maximum  crack  depth,  cm„r 
Default  =  thickness  of  plate. 

Note  -  This  analysis  does  not  require  a  BETA  solution  option. 

-  This  analysis  is  three  dimensional,  with  a  along  the  surface  and  c 
through  the  depth. 

b)  RETARD  option 

Line  1:  RETARD 

Line  2:  Cl  ,  C2  ,  C3  ,  C4  ,  C5  ,  CC 

Cl  -  defines  the  retardation  model  to  be  used: 

Cl  =  1;  Wheeler  retardation  model 
Cl  =  2;  Willenborg  retardation  model 

C2  -  indicates  to  the  program  whether  plane  stress  or  plane  strain  is 
to  be  assumed. 

C2  =  0:  plane  stress 
C2  =  1:  plane  strain 

Default  =  plane  stress 

C3  -  If  Cl  =  l  then  C3  is  Wheeler’s  retardation  parameter  m. 

-  If  Cl=2  then  C3  (if  not  zero),  represents  the  maximum  overload 
in  the  Gallagher-modified  Willenborg  model. 

C4  -  flags  whether  a  re-run  of  the  program,  without  retardation,  is  to 
occur. 

C4  =  0:  job  re-run  without  retardation 
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C4  =  1:  no  re-run. 

Default  =  job  re-run 

C5  -  gives  any  initial  overload  condition  of  the  specimen  at  program 
start. 

Default  =  0.0 

C6  -  gives  the  initial  plastic  zone  size  of  the  specimen  at  program  start 
due  to  a  previous  overload,  ap  . 

Default  =  0.0 

Note:  if  Ci=2  and  C3  ^  0  then  only  the  FORMAN  and  FORMANMOD 
equations  may  be  used. 

c)  BETA  option 

In  the  BETA  options,  the  upper  and  lower  crack  size  limits  define  the  range 
over  which  a  particular  P  solution  will  be  applied.  Thus  it  is  possible  to  have 
several  P  solutions  being  used  over  different  phases  of  the  crack’s  growth. 

Note:  if  some  of  the  values  (Cl,  C2,  ...  CO)  required  below  are  zero,  then 
they  must  be  entered  as  such.  Leaving  blanks  is  not  acceptable  and 
will  cause  the  program  to  crash. 

Line  1:  BETA 

Line  2:  Cl,  C2,  C3,  C4,  C5,  CO 

Cl  =  1:  Constant  correction  factor  model. 

C2  -  correction  factor 
C3  -  lower  crack  limit 

Default  =  initial  crack  size 
C4  -  upper  crack  limit 
Default  =  infinity 
C5  -  zero 
C6  -  zero 

Cl  =  2:  finite  width  correction  model 
C2  -  effective  plate  half  width 
C3  -  lower  crack  limit 

Default  =  initial  crack  size 
C4  -  upper  crack  limit 
Default  =  infinity 
C5  -  zero 
C6  -  zero 

Cl  =  3:  Tabulated  function  of  a/L  versus  p 
C2  -  half  width  of  specimen 
C3  -  n  -  number  of  points  in  table 
C4  -  lower  crack  limit 

Default  =  initial  crack  size 
C5  -  upper  crack  limit 
Default  =  infinity 
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C6  -  zero 

plus  Line  3:  Table  heading  (80  characters  maximum) 

Line  4:  l'1  point,  (a/L) i,  /?| 

Line  5:  2'“*  point,  (a/L) 2,  /?2 

Line  n  +  3:  n"‘  point,  ( a/L)„ ,  /?„ 

L  is  a  characteristic  length  which  for  CRKGRW  is  the  half  width 
of  the  specimen. 

Hence,  the  a/L  versus  /?  table  is  actually  a/(.(  width)  versus  ft 

Cl  =  4:  Exactly  the  same  as  Cl  =3.  The  reason  for  the  inclusion  of  two 
table  options  is  that  the  two  tables  can  be  applied  over  differ¬ 
ent  stages  of  the  crack’s  growth.  This  is  achieved  by  specifying 
different  upper  and  lower  crack  limits  (C4  and  C5). 

Cl  =  5:  Uniaxial  Bowie  solution  for  single  crack  from  a  circular  hole. 

02  -  radius  of  hole,  r. 

C3  -  lower  crack  limit 

Default  =  initial  crack  size 
C4  -  upper  crack  limit 
Default  =  infinity 
C5  -  zero 
C6  -  zero 

Cl  =  6:  Uniaxial  Bowie  solution  for  two  cracks  from  a  circular  hole. 

C2  -  radius  of  hole,  r. 

C3  -  lower  crack  limit 

Default  =  initial  crack  size 
C4  -  upper  crack  limit 
Default  =  infinity 
C5  -  zero 
C6  -  zero 

Cl  =  7:  ASTM  Compact  Tension  specimen 
C2  -  width  of  specimen,  b. 

C3  -  thickness,  t. 

C4  -  lower  crack  limit 

Default  =  initial  crack  size 
C5  -  upper  crack  limit 
Default  =  infinity 
C6  -  zero 

Cl  =  8:  Grumman  Compact  Tension  specimen 
C2  -  width  of  specimen,  b. 

C3  -  thickness,  t. 

C4  -  lower  crack  limit 

Default  =  initial  crack  size 
C5  -  upper  crack  limit 
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Default  =  infinity 
CO  -  zero 

Cl  =  9:  Lui  solution  for  a  double  quarter  crack  from  a  circular  hole. 

C2  -  radius  of  hole,  r. 

C-3  -  initial  crack  size,  atl . 

C4  -  final  crack  size 
C5  -  zero 
CG  -  zero 

Cl  =  10:  Newman  and  Raju  solution  for  a  corner  crack  extending  from  a 
hole.  Three  dimensional. 

C2  -  radius  of  hole,  r. 

C3  -  single  or  double  crack 
C3  =r  1 :  single  crack 
C3  =  2:  double  crack 
C4  -  initial  crack  depth,  c(). 

Default  =  thickness  of  specimen,  t. 

Co  -  half  width  of  specimen,  6. 

C6  -  thickness  of  specimen,  t. 

Cl  =  11:  Newman  and  Raju  solution  for  a  surface  crack  extending  from 
inside  circular  hole. 

C2  -  radius  of  hole,  r. 

C3  -  single  or  double  crack 
C3  =  1:  single  crack 
C3  =  2:  double  crack 
C4  -  initial  crack  depth,  Co. 

Default  =  half  thickness  of  specimen,  t. 

C5  -  half  width  of  specimen,  b. 

CG  -  half  thickness  of  specimen,  t. 

d)  END 

Line  1:  END 

This  keyword  is  required  to  terminate  an  ANALYSIS  segment.  It  must 
always  appear  in  the  data  file  after  all  the  analysis  options  have  been 
completed. 
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4.1.8  LOADS  segment 

This  segment  defines  the  load  sequences  for  one  or  more  missions.  Each  mission  is 
defined  as  a  sequence  of  loads  (or  stresses)  and  the  number  of  applications  of  each 
load.  The  order  in  which  the  various  missions  are  specified  is  not  important  (the 
order  in  which  they  will  be  applied  is  controlled  by  the  SPECTRUM  segment). 
The  program  numbers  them  consecutively  as  they  appear  in  this  segment. 

Line  1:  LOADS 
Line  2:  TITLE 
Line  3:  NBLKS  ,  LPRT 
Line  4:  MULT 

TITLE  -  comment  line  describing  load  sequence  (80  characters  maximum). 
NBLKS  -  maximum  number  of  times  the  whole  program  block  is  to  be  repeated 
(must  be  an  integer). 

LPRT  -  print  mission  specifications  to  output  file.(must  be  an  integer) 

LPRT  =  0:  print  information 
LPRT  =  1:  do  not  print  information 
MULT  -  multiplication  factor  on  all  loads  read  in. 

Note:  This  option  is  very  useful  as  it  allows  non-dimensional  spectra  or  strain 
gauge  spectra  to  be  used  and  converted  to  stresses  or  loads  with  the 
MULT  factor. 

There  are  three  ways  of  entering  the  loads  into  the  program.  Each  mission  starts 
with  a  keyword/header  line.  The  keywords  are:- 

MAX-MIN 

R-DELTA 

MEAN-ALT 

The  keywords  are  in  columns  1-10,  followed  by  an  optional  20  character  heading 
describing  the  mission  in  columns  11-30 

Lines  of  data  follow,  each  representing  a  single  layer  of  the  mission.  A  layer  consists 
of  a  pair  of  stress  parameters  followed  by  the  number  of  cycles  for  which  they  are 
to  be  applied. 

Line  5:  MAX-MIN 

Line  6:  1st  data  line,  0maT,  ,  <T„„n,  >  No.  of  cycles 
Line  7:  2n<1  data  line,  crn,alJ  ,  ,  No.  of  cycles 

Line  n-t-5:  n,h  data  line,  cr,„arn  ,  aw,„n  ,  No.  of  cycles 
Line  n+6:  -9999  -9999  -9999 

Note:  Line  n+6  is  used  as  a  termination  line  due  to  the  free  format  read. 
Line  5:  R-DELTA 

Line  6:  1st  data  line,  Aai  ,  R\  ,  No.  of  cycles 
Line  7:  2n<1  data  line,  Aa?  ,  R2  ,  No.  of  cycles 
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Line  n  +  5:  n'1'  data  line,  Ac,,  ,  R„  ,  No.  of  cycles 
Line  n+ 6:  -9999  -9999  -9999 

where  =  Ao/(l  -  R) 

and  <7 ,,,,,,  =  c,„rtJ.  Ac 

Line  5:  MEAN-ALT 

Line  C:  1st  data  line,  ,  o„;(l  ,  No.  of  cycles 

Line  7:  2ni1  data  line,  c„„vlBj  ,  c,,/^  ,  No.  of  cycles 

Line  n  +  5:  nth  data  line,  ,  ca/(ii  ,  No.  of  cycles 

Line  n+6:  -9999  -9999  -9999 

where  =  o ,,, „  +  c„p 

and  C,B/B  in  ran  &all 

Note:  When  the  BETA  options  7  and  8  are  chosen  in  the  analysis  section, 
the  layers  must  be  presented  as  loads  and  NOT  stresses. 

To  terminate  the  LOADS  segment  the  final  line  must  have  the  following  keyword:- 
END  LOADS 

Note:  At  present  CRKGRW  truncates  all  negative  loads  to  zero. 

4,1.9  SPECTRUM  segment 

This  segment  defines  the  pattern  in  which  missions  are  selected  from  the  data 
presented  in  the  LOADS  segment  to  make  a  complete  block. 

Line  1:  SPECTRUM 
Line  2:  NCYC,  HR/BL 
Line  3:  NSEGS,  IPRT 
Line  4:  NFLTSi,  MISS^ 

Line  5:  NFLTS2,  M1SSN2 

Line  NSEGS  +  3:  NFLTSxsecs>  MISSN^sehs 

NCYC  -  this  is  the  total  number  of  cycles  in  a  block 
HR/BL  -  this  is  the  total  number  of  hours  per  block. 

NSEGS  -  the  number  of  segments  in  the  block.  That  is,  the  number  of  lines  in 
the  spectrum  definition  to  follow  (must  be  an  integer). 

IPRT  -  print  mission  specifications  to  output  file  (must  be  an  integer). 

IPRT  =  0:  print  information. 

IPRT  =  1:  do  not  print  information. 

NFLTS  -  number  of  times  a  mission  is  to  be  repeated  (must  be  an  integer). 
MISSN  -  the  mission  number  from  LOADS  segment  (must  be  an  integer). 
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4.1.10  PRINT  segment  (optional) 

This  segment  controls  the  print  of  information  to  the  output  file. 


Line  I:  PRINT 
Line  2:  I  ,  J  ,  K  ,  L,  A 

I  -  crack  growth  output  every  I’th  block  (integer). 

J  -  crack  growth  output  every  J’th  segment  (integer). 

K  -  crack  growth  output  every  K’th  flight  (integer). 

L  -  layer-by-layer  crack  growth  output  every  L’th  flight  (integer). 

A  -  crack  growth  output  whenever  Aa  and/or  Ac  >  A. 

A  zero  corresponds  to  no  output  and  is  the  default. 

4.1.11  PLOT  segment  (optional) 

This  segment  allows  the  user  to  obtain  output  files  containing  lists  of  two  variables. 
The  output  files  will  have  a  ‘.PG[number]’  extension.  The  format  of  this  segment 
is  shown  below:- 

Line  1:  PLOT 
Line  2:  x-variable 
Line  3:  y-variable 

The  possible  keywords  for  the  variable  names  are  listed  below:- 


KEYWORD 

MEANING 

A 

a 

C 

c 

CYCLES 

HOURS 

DA/DN 

da/dN 

DC/DN 

dc/dN 

DELKA 

A K  in  o  direction 

DELKC 

A K  in  c  direction 

BLOCKS 

The  ‘.PG|number|’  file  will  contain  two  columns  of  numbers.  The  first  column  will 
bo  the  x-variable  and  the  second  column  will  be  the  y-variable.  These  files  can 
then  be  used  in  conjunction  with  inhouse  plot  routines.  The  interval  of  writes  to 
this  file  is  linked  to  that  of  the  PRINT  segment.  Care  must  be  taken  to  ensure 
that  huge  files  are  not  created. 

4.1.12  END  DATA  segment 

This  keyword  signals  to  the  program  that  the  input  has  ended. 


4.2  LIMITATIONS 

CRKGRW  has  the  following  limits:- 

-  Maximum  of  50000  layers  per  mission. 

-  Maximum  of  50  missions  per  block. 

-  Maximum  of  200  segments  per  block. 

-  Maximum  of  100  values  in  the  A K  ,  da/dn  table  (DA/DN  option). 

-  For  the  R-DA/DN  option,  the  total  number  of  lines  in  all  the  files  must  not 
exceed  100. 

-  The  name  of  the  data  input  file  must  not  exceed  24  characters  (including 
the  .DAT  extension  -  see  section  5). 

-  The  filenames  requested  by  the  DA/DN  and  R-DA/DN  options  must  not 
exceed  10  characters. 

4.3  COMMON  INPUT  FILE  ERRORS 

There  are  a  number  of  errors  which  may  occur  when  C'RKGRVV  roads  the  input 
data  file.  Two  common  errors,  along  with  their  causes  are  given  below. 

Error  CRKGRW  error  ‘Unrecognised  keyword’ 

Causes  -  Illegal  keyword 

-  Comment  line  without  a  ‘C’  at  the  beginning 

-  Using  lowercase  letters  in  the  keyword 
eg;  Surface,  Lines  Per  Page,  BeTA 
Mistyped  keyword 

eg:  PRN1T,  WALKER,  LINESPER  PAGE 

The  keyword  does  not  begin  in  the  leftmost  column  of  the  line  (column 

1) 

eg;  C  This  is  a  comment  line  with  the  ‘C’  in  column  1 

ANALYSIS 

C  The  previous  line  will  not  be  recognised  because  the 
C  ANALYSIS  keyword  does  iu>t  begin  in  column  1. 

Error  -  Program  crashes  while  trying  to  road  the  data  file 
Causes  -  Placing  a  real  number  where  an  integer  is  expected. 

eg;  LPRT  must  be  either  0  or  1,  however,  putting  0.,  (Ml,  1.  or  1.0 
in  the  input  file  will  cause  the  program  to  crash. 

-  not  putting  in  enough  numbers  as  required  by  the  option 

eg:  The  BETA  option  requires  G  numbers  to  be  on  the  line  after 
BETA,  but  some  of  the  sub-options  require  only  four  or  five  non¬ 
zero  values.  For  example,  if  Cl  -  2  is  chosen  then  only  C2,  CTj  and 
C4  are  non  zero  but  C5  and  C’G  must  still  be  entered  as  zeros. 

-  There  was  no  909!),  9999,  9999  line  at  the  end  of  a  set  of  layers  in 

the  LOADS  segment.  This  will  i.Mi.-v  tin  progiam  to  attempt  to  read 
the  next  keyword  as  a  number  which  causes  the  crash. 
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5.  RUNNING  CRKGRW 


The  input  file  which  contains  all  the  information  for  CRKGRW  must  have  a  .DAT 
extension  to  its  name  (eg:  F18LUG.DAT,  SPAR.DAT,  TEST.DAT).  The  name 
may  be  any  legal  combination  of  upper  and  lower  case  characters.  When  CRK¬ 
GRW  is  run,  it  will  ask  for  the  name  of  this  file.  The  name  should  then  be  entered 
in  uppercase  only.  You  do  not  need  to  type  the  .DAT  extension  since  CRKGRW 
automatically  searches  for  .DAT  files.  For  example,  entering  F18LUG  will  cause 
CRKGRW  to  read  data  from  F18LUG.DAT. 

CRKGRW  produces  several  output  files.  If  your  input  file  was  F18LUG.DAT 
then  the  output  files  produced  would  be 

FILE  FILE  INFORMATION 

F18LUG.EXT  indicates  if  extrapolation  occurs. 

F18LUG.OUT  main  output  file. 

F18LUG.PG1  1"  plot  output  file. 

F18LUG.PG2  2“a  plot  output  file,  etc 

The  program  requires  the  name  of  the  data  file  at  runtime,  which  it  requests. 
MATERIAL  options  DA/DN  or  R-DA/DN  may  require  separate  data  files.  These 
should  be  created  before  running  CRKGRW  and  their  names  entered  when  CRK¬ 
GRW  requests  them.  To  end  the  entry  of  MATERIAL  data  files  to  CRKGRW, 
the  keyword  ‘NOMORE’  should  be  typed  instead  of  a  filename. 

When  reading  through  the  data  input  file,  CRKGRW  will  stop  at  the  first  error 
it  encounters.  When  this  happens  you  must  edit  the  data  file  to  remove  the  error 
and  then  re-run  CRKGRW. 

An  example  data  file  is  presented  in  Appendix  1.  This  data  file  is  only  intended 
to  be  an  example  of  a  variety  of  options  which  were  presented  in  Section  4. 
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FIGURE  4:  DEFINITION  OF  a  AND  c  FOR  A  SURFACE  FLAW 


IGURE  5:  DEFINITION  OF  a  AND  c  FOR  A  SINGLE  OR  DOUBLE  QUARTER 
ELLIPTICAL  CORNER  CRACK  FROM  A  HOLE 


Dimension 

ASTM 

GRUMMAN 

in) 

2.2,  2.5 

2.5 

/b 

0.60 

0.95 

d 

/b 

0.25 

0.60 

FIGURE  7:  ASTM  AND  GRUMMAN  COMPACT  TEST  SPECIMEN 


Mission  1  -  defined  by 

stating  each  layer's 

o  ,  o  .  and  number 
max  min 

of  cycles 


Mission  2  -  defined  by 

stating  each  layer's 

o  ,  a  , .  and  number 
mean  alt 

of  cycles 


Mission  3  -  defined  by 

stating  each  layer's 

a  ,  o  and  number 
max  mean 

of  cycles 


Mission  4  -  defined  by 
stating  each  layer's 
R,  ;\u  and  number  of  cycles 


FIGURE  8(a):  LAYERS  AND  MISSIONS 
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FIGURE  8(b):  THE  MISSION  OF  FIGURE  8(a)  APPLIED  TO 
FORM  A  BLOCK 


APPENDIX  I 


This  appendix  gives  an  example  CRKGRW  data  file,  with  extra  options  as  exam¬ 
ples  of  different  inputs. 

C 

C  THIS  IS  A  COMMENT  LINE  EXAMPLE 

C 

TITLE 

•1 

THIS  IS  A  SAMPLE  PROGRAM  TO  DEMONSTRATE  THE  INPUT 
OF  INFORMATION  INTO  CRKGRW. 


EQUATION 

WALKER 

MATERIAL 

ALUMINIUM  ALLOY  7050  -  material  comment  line 

1.301E-10,  .55,  3.25 

50.0  03.0 

THRESHOLD 

3.0  .1 

LIMITS 

.01  2.0  0.0  0.0 
ANALYSIS 

C - using  a  two  crack  Bowie  solution 

BETA 

5  0.25  0.0  0.0  0.0  0.0 

C - with  the  finite  width  correction  factor 

BETA 

2  2.0  0.0  0.0  0.0  0.0 

END 

LOADS 

SPECTRUM  FULL  DATA 
100  1 
30.0 

MAX-MJN 


1.0 

-0.0 

10 

) 

1.2 

-0.7 

28 

) 

0.9 

-0.2 

19 

)  Mission 

0.3 

0.0 

100 

) 

0.5 

0.3 

2  GO 

) 

-9999 

-9999 

-9999 

MEAN- ALT 

1.0 

0.2 

300 

) 

o.s 

0.1 

300 

)  Minion 

0.8 

0.3 

300 

i 

-9990  -9909  -9990 


END  LOADS 
SPECTRUM 
11853  1000 

3  0 

7  i  )  Segment  1:  7  (lights  of  mission  1  ) 

9  2  )  Segment  2:  9  flights  of  mission  2  )  1  block 

2  1  )  Segment  3:  2  flights  of  mission  1  ) 

C - - print  every  layer 

PRINT 

0  0  0  1  0 

q - or  print  every  five  blocks  and  whenever  crack 

C- - growth  exceeds  0.05 

PRINT 

5  0  0  0  0.05 

PLOT 

A 

HOURS 
PLOT 
DA/DN 
HOURS 
END  DATA 
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